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Abstract 

Scenarios of large extra dimensions have enhanced the importance for the study of black holes 
in higher dimensions. In this paper, we analyze an axisymmetric system of two black holes. 
Specifically the Bowen-York method is generalized for higher dimensions in order to calculate the 
initial data for head-on collision of two equal-mass black holes. Then, the initial data are evolved 
adopting the close-slow approximation to study gravitational waves emitted during the collision. 
We derive an empirical formula for radiation efficiency which depends weakly on the dimensionality. 
Possible implications of our results for the black hole formation in particle colliders are discussed. 

PACS numbers: 04.70.-s, 04.30.Db, 04.50.+h, ll.25.-w 
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I. INTRODUCTION 

Motivated by scenarios of large extra dimensions 

HQS 

, the study of black holes in 

higher dimensions has attracted attention. In these scenarios, extradimensional effects play 
an important role for the properties of black holes with radius smaller than the size of extra 

lanned acceler- 
ators. If the Planck energy is of O(TeV) as suggested in scenarios of 



dimensions. An interesting phenomena is mini-black hole production inpiarj 

m 

associated with quantum gravity will be observed at the CERN Large Hadron Collider (LHC) 



phenomena 



through black hole production 



HE 



5|. At the LHC, the produced black hole is expected to 
settle down to a higher- dimensional Kerr black hole after emission of gravitational waves. 
The formed black hole will be subsequently evaporated by the Hawking radiation, for which 
high-energy particles are emitted and may be detected. To know the parameters of the 
formed Kerr black hole which determines the feature of the Hawking radiation observed, 
studies for gravitational radiation are an important issue. 

In this paper, we study axisymmetric collision of two black holes in higher dimensions as 
a first step for the investigations on the black holes formed in particle colliders 1 . This is 

n 

an extension of our previous study [6( of the higher-dimensional Brill-Lindquist initial data 
and its temporal evolution by the close-limit approximation. In four dimensions, there is 
the well-known Bowen-York method [?| for generating initial data of several moving black 
holes. We generalize this method for higher dimensions. Similarly to the four- dimensional 
case, the conformally transformed extrinsic curvature is given analytically and the conformal 
factor should be calculated numerically. We perform numerical calculation for the conformal 
factor in the case of the head-on collision of two equal-mass black holes and then determine 
the location of apparent horizon (AH) that encloses two black holes (the common AH). We 
clarify the parameter space for common AH formation. 

Then, the initial data is evolved adopting the close-slow approximation. In this approxi- 
mation, the distance between two black holes and the linear momentum of each black hole 
are regarded as small parameters compared to the black hole size and the gravitational mass 
of the system, respectively. In this case, the system may be considered to be a perturbed 



1 Here, the incoming black holes are regarded as substitutes for the incoming particles. We do not consider 
the collision of produced black holes, since they are evaporated within the time scale of 10 -27 s but only 
one black hole is produced per 1 s. 
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single Schwarzschild black hole. In the four-dimensional case, the results in the close-slow 
approximation agree with those in the fully nonlinear analysis (see js, [9] and references in 
6j), and hence, it is natural to expect that this is also the case in the higher-dimensional 
case. We provide a formula for the radiation efficiency (i.e., the ratio of the radiated energy 
to the system energy) and discuss their dependence on the dimensionality. 

The purpose of this paper is twofold. One is to present a formulation for providing the ini- 
tial condition of colliding black holes in higher dimensions. Gravitational waves at the colli- 
sion of two black holes are accurately computed only by a fully nonlinear numerical relativity 
simulation. Although such a simulation has not been done yet for higher- dimensional space- 
time, several groups have developed robust techniques for numerically computing merger of 
binary black holes in four dimensions [lOj]. The methods used in four dimensions can be 
applied for higher- dimensional problems and hence, the simulation will be done in the near 
future. For such future simulations, a method for providing initial condition we present here 
will be useful. 

The second purpose is to approximately evaluate gravitational radiation in the collision 
of two black holes in the higher- dimensional case. Since we adopt a linear perturbative 
method, a reliable result is derived only for a small parameter space. However, our result 
for such parameters will be useful for calibrating fully nonlinear results in the near future. 

This paper is organized as follows. In the next section, we derive a generalized Bowen- 
York formulation |7| for higher-dimensional space. Applying the formulation, we construct 
the initial data for head-on collision of two equal-mass black holes in Sec. III. The ADM mass 
and the common AH are analyzed. In Sec. IV, we evolve the initial data adopting the close- 
slow approximation. The master variable of the linear perturbation around a D-dimensional 
Schwarzschild black hole is calculated numerically and the formula for radiation efficiency is 
derived. The radiation efficiency is shown to depend weakly on the dimensionality. Sec. V 
is devoted to a discussion on possible implication of our results for the radiation efficiency in 
the particle collision. In Appendix A, we present some solutions of the extrinsic curvature 
in the generalized Bowen-York formulation that were not introduced in Sec. II. In Appendix 
B, the method for imposing the initial condition for the master variable is explained. 
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II. HIGHER-DIMENSIONAL BOWEN-YORK METHOD 



In this section we generalize the Bowen-York formulation for higher dimension space and 
present a method for generating initial condition of N black holes with linear momenta. 

A. Formulation 

Let S(/i^, K^) be a (D — l)-dimensional spacelike hypersurface £ with the metric 
and the extrinsic curvature K^ u in a D-dimensional spacetime. We introduce a number 
n = D — 2. The Hamiltonian and momentum constraints are 

(n+i) R _ k^K^ + K 2 = 0, (1) 

V^K^ - h^K) = 0, (2) 

where ( n+1 ">R denotes the Ricci scalar of E and V M is the covariant derivative with respect 
to hfa,. Introducing defined by 

V = ^ {n -^ V, (3) 

the Hamiltonian and momentum constraints are rewritten as 

n — 1 



An 



y /^2{n+\)/{n-l) K u\ _ ^2(n+l)/(n-l)y ^ ?_^(n+3)/(n-l) y ^ = g /gN 

^ V M/ ^ Tl 1 

where and V M denote the Ricci scalar and the covariant derivative with respect to 

h^, respectively. Here raising and lowering the index are done with h^. 

In the following, we assume the conformal flatness on £, h^ v = 5^, and impose the 
maximal slicing condition, K = 0. A weighted extrinsic curvature is defined by K^ u = 
and hereafter its index is raised and lowered by 5^ (i.e., Kg = ^r 2 ( n+1 )/( n - 1 ) anc [ 
K^ u = ip r2 ( n + 3 )/( n - 1 ) K^ u ). Then the Hamiltonian and momentum constraints become 

d^K^ = 0, (7) 
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where denotes the ordinary derivative with respect to the Cartesian coordinate and 
Vf = <9 M <9 M . Following Bowen and York we assume that K^ v does not have the tensor 
mode and thus takes the following form: 

K, u = d,W v + d v W^ - -^—S^dpW. (8) 

n + 1 

Substituting this formula into the momentum constraint (|7j), we obtain 

77 — 1 

V f 2 ^ + —-dpd v W = 0. (9) 

71 + 1 

Introducing auxiliary functions _B M and Xi we decompose as 

Stj 4- 1 

W, = —-B, - {d, X + x u d,B v ) . (10) 



Then the equation becomes 

= ^V f 2 B M - -^rdtfx - * V d^\B v - n -^-d, (x^B v ) . (11) 
n — 1 n + 1 n + 1 v 

Hence the momentum constraint is satisfied if 



V 2 5 M = 0, (12) 
V f 2 x = 0. (13) 



Since the solutions of Eqs. (112]) and (113]) are analytically given, solutions for K ab are easily 
provided. 



B. iV-black-hole solutions 



To give linear momenta of black holes, we choose the solution 



where P M is a constant vector, G the gravitational constant, fl n the n-dimensional area of a 
unit sphere, and R = Then, we obtain 

^ = ^o^!^ + p ^ + P p^ P [(« " - V]} > (15) 

where n M = x^/R. This solution provides the extrinsic curvature for one boosted black hole 
located at R = 0. Actually, P M agrees with the ADM momentum: 

P m = ^ / - ^V) ^ (16) 



Since the momentum constraint ((Tj) is a linear equation, we can superpose N solutions. 
Denoting the locations of iV black holes as x£(a = 1,...,N), a solution of the extrinsic 
curvature is written as 

N 

471 

~nQ n Rl 

(17) 

where R a = \x^ — x%\, = (x M — x£)/R a , and (P a )n denotes the momentum of the ath 
black hole. 

The conformal factor ^ is obtained by solving the Hamiltonian constraint (jH]). Following 



a=l 



11], we assume that \P has the following form 

& = &BL + i>, (18) 

where 

^bl-I + ^E]^! (19) 
and M a denotes the mass parameter of ath black hole. Then, the equation for ip becomes 

Vf 2 ^ = -l!^il^^(^ BL + ^)-(3n+l)/(n-D_ (2Q) 

The solution in this procedure represents the so-called "puncture" space with N Einstein- 
Rosen bridges and N + 1 asymptotically flat regions (say, one upper sheet and N lower 
sheets) . 

Since the right hand side of Eq. (1201) behaves like 0(i?" _1 ) for R a — > 0, there is a regular 



solution for ip, which can be solved numerically (this fact was first pointed out in 11[ for 
the four- dimensional case). The ADM mass Madm is given by 

M ADM = —, tt^ / d^n^dS. (21) 

4vr(n - 1)G J R ^ 

Using the Gauss law, we find 

N 

M ADM = Y M a + — — / K^K^-^ n+1 ^ n ^d n+1 x. (22) 



ti 16nG J* 



The first and second terms could be interpreted as the sum of the mass of N black holes 
and the kinetic energy of the black holes, respectively. 
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III. INITIAL DATA FOR HEAD-ON COLLISION OF TWO BLACK HOLES 



In this section, we present initial data of the axisymmetric two-black-hole system following 
the formalism described in the previous section. 



A. Calculating conformal factor 

We introduce the Cartesian coordinate (z,Xk), (k = 1, ...,n) for (n + l)-dimensional space 



and write the solution of K^ v as 



K IV = K^ + K^\ (23) 



1 flU ±s -fJ,U 1 fiU 

where two black holes (±) are located at (z,Xk) = (±zq,0) and have momenta 
(=pP, 0). We assume that the two black holes have the same mass M + = M_ = Mo/2. 
The gravitational radius of a black hole of mass M is defined by 



16nGM 



l/(n-l) 



r * (M °' = l^arJ ' (24) 

In the case z = 0, the common apparent horizon (AH) that encloses two black holes is 
located at R = Rh(M ) where 

R h (M ) = A- l ^r h (M ). (25) 



Using Rh(Mo) the conformal factor in the puncture framework [U| is given by 

W = V BL + V, (26) 

where 

^bl = 1 + \ [R h {M a )r l + -Lj) . (27) 

For a numerical solution of ip, we introduce the cylindrical coordinate (z, p) where p = 
vSIi=i x t- m this coordinate, the equation for ijj becomes 

^, PP + ^ + { -^—^-t P + ^K^k^-^^ = 0. (28) 
p An 

Since the system is axisymmetric and equatorial-plane symmetric, it is sufficient to solve 
the equation for the domain < p < p max and < z < z max with the boundary conditions 
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i/j :P = at the z-axis and ip jZ = at the equatorial plane. For it! — > oo, ip asymptotically 
behaves as 

4nG(M ADM - M ) 



+ 0(1/R n ). 



(29) 



nVt n R n - 1 

Thus, we impose the so-called Robin condition, ip t R = — (n — l)ip/R, at the outer boundary. 
In the cylindrical coordinate, it is rewritten as 



^,pP + ^j, z z = -(n-l)ip. 



(30) 



We put the outer boundary at p max = z max = 5i?^(M ) and solve the equation for ip using 
a finite difference method with the grid number (101 x 101). Numerical computation was 
carried out for the parameter space < z /R h (M ) < 1 and < P/M < 1 with the 0.1 
interval or the 0.01 interval. 

There are two sources for numerical error of the conformal factor: One is associated with 
the finiteness of the grid spacing and the other with the finiteness of the outer boundary 
location. To evaluate the numerical error by the grid spacing, we took reference data with 
sufficiently large grid number (401 x 401 grids) and evaluated following three characteristic 
numerical error values: 



ei = 

e 2 = max 
e 3 = max 



Sat 


1 i (ref) 









(ref) 



(ref) 
N 



(31) 



4'" 



)■ 



where iV stands for the label of the grids. Changing the grid spacing, we confirm that 
€i decreases with improving the grid resolution at second order. All the three error esti- 
mates were found to be small. With our standard choice, e 2 is 0.005-0.2% for D = 4- 
11. Computation was also performed changing the location of the outer boundaries from 
Pmax = Zmax = 5i?/i(M ) to 20Rh(M ) and the numerical solution converges with increasing 
the radius of the outer boundaries. In this case, e 3 is found to be relatively large, about 
9-17% for D = 4-11. Such large differences occur at outer boundaries and come from the 
fact that the Robin condition is an approximate boundary condition. However, because ip is 
small at the outer boundary, the error of the conformal factor W is small. In fact, we found 
that e 2 is about 0.3% for D = 4 and becomes smaller as D is increased (less than 0.01% in 
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the case D = 11). The error e\ is smaller than £3 by a factor for D = 4 and is comparable to 
62 for D = 11. Because ^ is used in the calculations of the ADM mass or the AH, the error 
for these calculations is expected to have the order of €2- Hence, we consider that sufficient 
accuracy is obtained in our calculation. 

Figured] shows the contours of Madm/M on the (zo/Rh(M ), P/M )-plane for D = 4-11. 
The difference between Madm and M indicates the strength of nonlinearity due to the right 
hand side of the Hamiltonian constraint ([6]). We find that the nonlinearity of the system 
becomes large as D is increased for fixed values of z Q /R h (M ) and P/M . 



B. Common apparent horizon 

A common AH that covers two black holes is searched for introducing the spherical 
coordinate (R,9) by R = a/ p 2 + z 2 and tan# = z/p and assuming that the location is 
denoted by R = h(6). Then the equation for the AH is given by 

V M s" -K + K^s v = 0, (32) 

where s M is the unit normal to the surface: 

(7/-2/(n-l) 

= , (1, -h e /R 2 ) . (33) 



1 + h 2 e /R 2 

Then, the equation for the AH reduces to an ordinary differential equation for h{9): 
2 Vr 



hao-n 



n-l & 



_ &-2n/(n-i) h ^ h2 + h 2 e (% RR _ 2K R ^ + Kee^j = 0. (34) 



This equation is solved under the boundary conditions h t $ — at 9 — and n/2. 

Figure [2] shows the critical line for the common AH formation. Irrespective of the dimen- 
sionality, the common AH is formed if Zo/Rh(Mo) is sufficiently small. It is interesting to 
point out that, for D = 4, the AH formation is more subject for larger values of P while 
for D > 5, the presence of P tends to prevent formation of the common AH. This might 
seem strange because the kinetic energy is naturally expected to help the AH formation. 
In order to obtain a corroboration for this result, we direct attention to the AH (denoted 
by R = h(6)) at the critical separation z = z^"^ in the case P = 0. Then we add small 




FIG. 1: The contour line for M A dm/M on the (z /R h (M ), P/M )-plane. 



10 



0.8 



0.6 



0.8 



0.6 



P/M 



0.4 



0.2 



0.8 



0.6 



P/M 



0.4 



0.2 



0.8 



11.6 



P/M 



0.4 



0.2 



»=4 



P/Mq AH formation 

0.4 

0.2 



0.2 0.4 0.6 0.8 1 
Zo/R h {M Q ) 
D=6 



No AH 



AH formation 



0.2 0.4 0.6 0.8 1 

zo/R,,(M ) 
0=8 



No AH 



AH formation 



0.2 0.4 0.6 0.8 1 

zo/R h (M ) 
D=10 



No AH 



AH formation 



0.2 0.4 0.6 0.8 1 

zo/RiAMo) 



O.N 



0.6 



P/M 



0.4 



0.2 



O.N 



0.6 



P/M 



0.4 



0.2 



O.N 



0.6 



P/M 



0.4 



0.2 



0.N 



0.6 



P/M 



0.4 



0.2 



D=5 



AH formation 



No AH 



0.2 0.4 0.6 0.8 1 

zo/ff/,(Mo) 
D=l 



No AH 



AH formation 



0.2 0.4 0.6 0.8 1 

zo/Ri,(M ) 
D=9 



No AH 



AH formation 



0.2 0.4 0.6 0.8 1 

Zo/Rl,(M ) 
D=ll 



No AH 



AH formation 



0.2 0.4 0.6 0.8 1 

zo/Ri,(M ) 



FIG. 2: The critical line for the common AH formation on the (zq / Rf l (Mo) , P/Mo)-plane. 
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FIG. 3: The expansion + of a minimal surface in the P <C Mq case at the critical separation 
zq = Zg Cnt ^. The cases D = 4 (left) and 6 (right) are shown. In the case D = 4, + is negative for 
all 9 and the minimal surface is a trapped surface. But in the case D = 6, there is a region of 
where #_|_ becomes positive, which means that the minimal surface is not a trapped surface. 

P <C Mq without changing the value of zq. Because the change in the conformal factor 
is 0(P 2 /Mq) and can be ignored, the surface R = h(9) is interpreted as a minimal sur- 
face on which V M s M = holds. Now let us consider the expansion on this minimal surface 
which is given by 9 + = K VbV s ii s v . Figure [3] shows the behavior of 9 + /P in the four- and 
six-dimensional cases. In the case D = 4, we find that 9 + is negative and thus the minimal 
surface is a trapped surface. This guarantees the existence of an AH outside of the minimal 
surface. Hence the motion helps AH formation in the four-dimensional case. On the other 
hand, in the case D = 6, 9 + becomes positive on some part of the minimal surface. The 
similar behavior was found for all 5 < D < 11. Therefore, the minimal surface is not a 
trapped surface and the motion does not necessarily help the AH formation in the higher- 
dimensional cases. We found numerically that the term KQgs 6 s e mainly contributes to the 
positivity of 9 + . This is because for high D, the AH at the critical separation is hourglass 
shaped in the neighborhood of the equatorial plane and the value of h t $ is quite large around 
9 ~ 1. It enhances the value of Kogs d s e which is proportional to h 2 e /h A . 
The AH mass Mah is determined by the AH area Aah as 



The area theorem of black hole constrains that the AH mass never decrease 2 . Thus, the 



2 Since the proof of the area theorem is not sensitive to the spacetime dimension, it holds also in the current 




(35) 
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TABLE I: The values of 1 - M All /M ADM evaluated on the AH critical line for P/M Q = 0.0, 0.5, 1.0. 
The unit is %. 
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7.7 
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7.4 


6.9 


P/Mq = 0.5 


2.8 


9.0 


12 


14 


16 


17 


18 


19 


P/Mq = 1.0 


5.5 


15 


19 


21 


22 


24 


25 


26 



AH mass provides the lower bound on the black hole mass of the final state. In other words, 
the quantity Madm — ^ah is the upper bound on gravitational radiation energy in the 
collision process. Table [I] shows the value of 1 — Mah/^adm on the AH critical line for 
P/M = 0,0.5, 1.0. As the value of D increases, 1 — Mah/^adm becomes large for fixed 
values of P/Mq = 0.5 and 1.0 on the AH critical line. The area theorem provides a stricter 
condition for smaller value of D. 



IV. CLOSE-SLOW ANALYSIS 



Next we compute gravitational waves in the head-on collision of two black holes adopting 
the close-slow approximation. In this approximation, we assume that Zq <C r/j(M ), P <C M , 
and zo/rh{M ) ~ P/M , and evaluate gravitational wave energy up to order of (z /rh(Mo)) 2 
using a linear perturbative approach. In the following, the gravitational radius of the system 
r h (M ) is used as the unit of the length [i.e., Th{M ) = 1] unless specified. 



A. Close-slow form of the initial data 

Since we analyze gravitational waves in the Regge- Wheeler type method, Kffl in the 
spherical-polar coordinate (R, 9, 0*), = fa, n -i)> should be derived: 



4M ~ R{±) T 2 

-K R = — —r COSU{R =F Zq cost 



(n + l)P H Rl +l 



+ Z ° T Z7 9 [(" -WT Zo cos 6f - Rl] , (36) 



system as long as naked singularities do not exist. 
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4M ±1 
(n + l)P R R n+1 



sin 9{R =F 2zq cosi 



and 



+ Z ° T J! + T 6 [±fa - l)*o sm 9(R ? z cos g)] , (37) 

4M - e(± ) 2 . 2 Z ^P C OS# r 2 . 2 21 mm 

(n + l)P e = RY 1 * Sm ^ ^ ~ ° Sm ±J ' ( > 

4M ^ t (±) _ zq T R cos g 



{n + l)P ^ R 
For z Cl, K^ v = Kffl + is expanded as 

Kg = (z P/M )^± [ n - 2 -(n 2 + n-2) cos 2 9] R~^ + O(z 2 P/M ), (40) 

Kl = ( Zo P/ Mo )^ [(n - 1) cos 2 + 1] + O(z 2 P/M ), (41) 

= (z P/M Q )^±I [(n + 1) cos 2 (9 - 1] P~( n+1 > + O(z 2 P/M ), (42) 

and i^Rg = O(2qP/M ). The leading-order term of K ab is found to be O(z Q P/M ) and 
hence the right hand side of the Hamiltonian constraint ([6]) is of order 0(zqP 2 /Mq). In the 
close-slow approximation adopted here, such terms are higher order and we ignore them. 
Thus, if) = and Madm = Mq in this approximation. 

As a result, the conformal factor is given by the Brill-Lindquist one: 

^ BL = l + i(^ + ^). (43) 

By transforming from the isotropic coordinate to the Schwarzschild-like coordinate 

r = R^ n ~ 1 \ ^0 = 1 + ^1, (44) 



we find that the system is regarded as a perturbed Schwarzschild black hole 
ds 2 - 





4/( 


% J 




\ 4/(n 


-1) 



(/ ' J + r 2 (^ 2 + sin 2 M 2 „ 1 ) 



= 1 + VC"-/)^ 1 m 2 cr^%os9) + 0(4), (46) 
!Pb / 1 + 1/4P™- 1 VP/ 2 v ; v o;> v ; 

where C*| A ' denotes the Gegenbauer polynomials defined by the generating function 

oo 

(l-2xt + t 2 )- x = Y^C l e X \x)t i . (47) 

1=0 
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Note that the metric (HST) with (|46|) is the same as that in our previous analysis of the time- 
symmetric initial data [6]. However, the time-asymmetry is present because of the presence 
of nonzero K^ v . 

As found above, the order of the perturbation of the initial metric is O(zq) and of the 
extrinsic curvature is O(z P/M ). In the following, we consider the situation where both 
zq and P/Mq have the same order. Under this condition, we can evolve the system using a 
standard perturbation method in the Schwarzschild spacetime. From Eqs. fj40l) - fj42l) . fj45l) . 
and (T46|) . the leading order of the perturbation contains only the £ = 2 mode. 



B. Time evolution by the master equation 

The gauge-invariant method for the perturbation around the Schwarzschild black hole 

□ 

was developed by Kodama and Ishibashi [12J . They derived a master equation for a variable 
which is related to the gauge-invariant quantities of the perturbation, as 

<9 2 $ <9 2 $ 

V s $ = 0, (48) 



dt 2 drl 
where 

Vs (r) = f{T)Q{r) (49) 
S[) 16r 2 tf 2 (r)' 1 ' 

and 

H{r) =m + {l/2)n(n + l)x, x = l/r n ~\ (50) 
m = k 2 - n , k 2 = £(£ + n- 1), (51) 

Q(r) = n 4 (n + l) 2 x 3 + n(n + 1) [A(2n 2 - 3n + 4)m + n(n - 2)(n - 4)(n + 1)] x 2 

- Yin [(n - 4)m + n(n + l)(n - 2)] mx + 16m 3 + An{n + 2)m 2 . (52) 

r* denotes the tortoise coordinate defined by 

m 

Initial values of $ and $ (a dot denotes the time derivative hereafter) are related to 
the metric perturbation and K ab , respectively. We describe the detail in Appendix B. The 
equation for $(0,r) is the same as that in the Brill-Lindquist case {&] and the solution is 
given by 

* (0 r) ( z 2 ) n v^K + 3n + 4 + ^ + 3)v7] 
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FIG. 4: Time derivative of the master variable $ in the unit of ZqP/Mq for D = n + 2 = 4, 6, 8, 10 
(left) and for D = 5,7, 9, 11 (right). For P > 0, <3? is negative. 

where the definition of is given in Eq. (1B3I) . This has the order of £g. O n the other hand, 
$(0, r) is proportional to ZqP/Mq and its value is obtained by solving Eq. (IB27I) described 
in Appendix B. The solution is 

2n ^/J 2{n + 2) + (n + l)x 



<S>(0,r) = -(z P/M ) 



(n - l)K [ 2 ] r n/2+1 2(n + 2) + n(n + 



(55) 



We show the behavior of $(0, r#) in Fig. HI 

Since Eq. f )48|) is linear, $ is naturally decomposed into two parts 

$ = (^)$BL + (ZqP/Mo) $BY- 



(56) 



The solution for <3>bl is t ne same as that in the time-symmetric case derived in our previous 
jg] 3 . Here, we show only the computation of $by- For_the numerical computation, 



6[] and solve the equation in 



paper 

we use the second-order finite differencing code developed in 
the domain —200 < r* < 1000 with the grid spacings dr* = 0.01 (D = 4-7) and 0.005 
(D = 8-11) and dt = 0.2dr*. Computation was performed changing the grid spacing and 
we confirmed that the numerical results converge at second order. For the chosen grid 
spacing, the error evaluated with e x in Eq. (J3TJ) is 0.05%-0.9% for D = 4-7 and 0.4%-2% 
for D = 8-11. The error in the values listed in Table II is < 0.1%. 



3 In [J $(0, r*) was solved numerically. We recalculated the temporal evolution using the analytic formula 
(l54l) of the initial condition and found that the difference between the two is ~ 10 _5 %. 
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FIG. 5: The time evolution of <£>by for D = 4-11 observed at r* = 100. 
C. Numerical results 



Figure shows numerical results for the time evolution of $by- Soon after the onset of the 
calculation, a quasi-normal mode is excited irrespective of dimension and, subsequently, the 
power-law tail is seen for four and odd dimensions. We read off the quasinormal frequencies 
cjq N and checked the consistency with previous results of c«Jq N listed in {(], Q, 14 1 . The 
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FIG. 6: The time evolution of $by (dashed line) and $bl (solid line) for D = 6 (left) and 10 
(right) observed at r* = 100. 

TABLE II: The values of a, c 2 and c 3 of Eq. (EH) for D = 4-11. 



D 


4 


5 


6 


7 


8 


9 


10 


11 


ci 


0.0252 


0.0245 


0.0290 


0.0288 


0.0258 


0.0223 


0.0194 


0.0172 


C2 


-0.165 


-0.243 


-0.294 


-0.287 


-0.251 


-0.213 


-0.182 


-0.158 


C3 


0.343 


0.671 


0.808 


0.765 


0.647 


0.539 


0.456 


0.396 



behaviors of $bl and $by ar e compared in Fig. [6] for D = 6 and 10. It is found that phases 
of $bl and $by disagree and the phase shift is ~ n. This implies that two terms interfere 
each other. This feature universally holds irrespective of the dimensionality. 

From $bl and ^by, we calculate the radiated energy of gravitational waves by the fol- 
lowing formula (see p, [ijj] for a derivation): 

k 2 {n-l){k 2 -n) 



rad 



327TOG 

Substituting Eq. ( 1561) into the above formula, E ra <i is rewritten as 

-Erad 



Cl zl + c 2 z 6 (P/M ) + c 3 z 2 (P/M y 



(57) 



(58) 



where Ci, c 2 , and C3 are constants determined by numerical integration. These values are 
listed in Table [TTJ The formula (158"]) together with Table [III will be used for the benchmark 
of the fully nonlinear analysis in numerical relativity as in the four-dimensional case 8j, |9| . 

Figure [7] shows the contours of i? ra d/M on the (zo/Rh(Mo), P/M )-plane. For a fixed 
value of z , the radiated energy decreases as the value of P/M is increased for small values 
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D=4 D=5 




zo/R h (M ) zo/R h (M ) 

FIG. 7: Contours of radiation efficiency E^/Mq (shown in the unit of %) on the 
(zo/Rf l (Mo),P/Mo)-p\ane for D = 4-11 predicted by the close-slow analysis. Note that the unit 
of zq is Rh(Mo) in this figure. 
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TABLE III: The values of E rad /M evaluated on the AH critical line for P/M = 0.0, 0.5, 1.0. The 
unit is %. 



D 


4 


5 


6 


7 


8 


9 


10 


11 


P/M = 


0.0 


0.0034 


0.059 


0.20 


0.34 


0.44 


0.49 


0.51 


0.52 


P/M = 


0.5 


0.30 


1.76 


2.83 


3.01 


2.80 


2.53 


2.28 


2.08 


P/M = 


1.0 


1.6 


7.3 


11.5 


12.2 


12.0 


11.4 


10.5 


9.8 



of P as P/M < 0.2 x zq. This is because the phases of $bl and $by disagree by a factor of 
« 7r (cf. Fig. 6). In the range P/M > 0.2 x zq, the amplitude of second term in Eq. ( 1561) 
exceeds the first term and _E ra d increases as P/M is increased. The similar behavior was 
reported in the four dimension case 9]. 



D. Dependence on dimensionality 

To get some insight for the dependence of radiation efficiency on the value of D, we 
evaluate E ra< x/M Q by choosing characteristic values of z and P for D = 4-11. In comparison, 
we fix the value of P/Mq since P/2Mq could be interpreted as the value of momentum divided 
by the rest mass, i.e., v/y/1 — v 2 where v is the velocity of each incoming black hole. We 
adopt P/Mq = 0,0.5 and 1 which corresponds to 7 = 1, v2 ~ 1.41, and \/E ~ 2.24. Then 
we recall Fig. [2] and evaluate -E ra d on the AH critical line, since the close-slow approximation 
holds for the system sufficiently close to the Schwarzschild spacetime and it is necessary to 
choose zq for which the common AH presents. 

Table [HI] shows the values of E Tajdl /MQ evaluated in this procedure. In the higher- 
dimensional cases, it is ~ 2%-3% for P/M = 0.5 and ~ 10% for P/M = 1.0. These 
values do not depend significantly on D in the higher-dimensional cases. Note that the 
values in Table II III does not contradict the bounds derived from the AH area shown in 
Table E 

The difference between Mo and Madm increases with increase of P/Mq (cf. Fig. [1]). This 
difference comes from time asymmetry of the initial data and has the magnitude of order 
(zqP/M ) 2 . Thus for a large difference with Madm/A^o — 1 ^ 0.1, the linear approximation 
breaks down. In the four- dimensional case, we know that the nonlinear effect suppresses 
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-E ra d/MADM [9|- If this is also the case for D > 5, the value of E ra ^/ 'Madm ~ 10% for P/M 
would be an overestimate and would indicate that gravitational radiation do not significantly 
contribute to the loss of the system energy in the head-on collision of two black holes. 



V. DISCUSSION 



In this paper, we presented a generalized Bowen-York formulation for the initial value 
problem of multi-black holes in higher- dimensional space. Using this formulation, we derived 
analytic solutions of the momentum constraint equation for extrinsic curvature of moving 
black holes in the conformally flat higher-dimensional space. As an application, initial data 
for head-on collision of two equal-mass black holes were computed numerically solving the 
Hamiltonian constraint for the conformal factor. The properties of the obtained space such 
as the ADM mass and the common AH were analyzed. We determined the critical line for 
formation of a common AH in the parameter space of (zq, P)- 

Then, we evolved this system using the close-slow approximation and analyzed gravita- 
tional waves. We derived the formula for the radiated energy for a wide variety of D [see 
Eq. (j58p together with Table HI] and clarified that it depends weakly on D; irrespective of 
D the efficiency is less than ~ 10% for P/M < 1 and Zq < r^(M ). 

As we mentioned in Sec. I, one of the motivations of this study is to clarify the total 
radiated energy of gravitational waves in the black hole formation through particle collisions 
in accelerators. However, there are two shortcomings in this study. One is that we use the 
linear perturbation method which provides a reliable result only for small values of z and 
P. The other is that the Bowen-York black holes system does not accurately model the 
high-energy particles system. 

As for the first problem, we can learn from the four- dimensional study. In , the results of 
the close-slow approximation and of the fully nonlinear numerical simulation are compared. 
It shows that the results in the close-slow approximation agrees approximately with the fully 
nonlinear results for small separation and small momentum. However, the difference becomes 
larger as separation or momentum is increased. For P ^> M , the nonlinear numerical 
results indicate that the radiation efficiency increases as P is increased and approaches to 
an asymptotic value, while extrapolation of the linear perturbation results gives infinity. To 
know the results for P ^> Mq in higher-dimensional spacetime, obviously, fully nonlinear 
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numerical simulation is necessary. 

The difference between the high-energy particles and the Bowen-York black holes is the 
other problem. As the model of a high-energy particle, the metric of Aichelburg-Sexl 15| is 
often used. It is obtained by boosting the Schwarzschild black hole to the speed of light with 
fixed energy p = 7m and taking lightlike limit 7 — > 00. The gravitational field is localized in 
the transverse plane to the direction of motion and forms a gravitational shock wave, which 
is a reminiscent of the infinite Lorentz contraction of the isotropic gravitational field of the 
original Schwarzschild black hole. The spacetime is flat except at the shock wave and this 
enables us to write down the metric of two high-energy particles outside of the lightcone of 
shock collision. Hence, the AH formation at the instant of collision has been investigated so 
far flHQ. 

Obviously, gravitational field (the conformal factor or the extrinsic curvature) is not 
localized in the case of Bowen-York black holes. This indicates that the Bowen-York black 
holes system cannot be a precise model for the high-energy particles system even in the case 
zq ^> Rh(Mo) and P Mq. Actually, there exist several discussions 19( about the fact that 
the Bowen-York black hole system contains unphysical gravitational waves, i.e., so-called 
junk energy. Thus, in the fully nonlinear numerical calculation with the Bowen-York initial 
data, one should take into account the possibility that junk energy changes the estimate of 
radiation efficiency To avoid this, one should make the junk energy radiate away before the 
collision setting the initial separation sufficiently large. 

Finally, we comment on the importance of the study of the grazing collision. In the head- 
on collision case, the upper bound on the radiation efficiency (< 40%) in the head-on collision 
of Aichelburg-Sexl particles was obtained [16| by studying the AH area and applying the 
area theorem. Although the same discussion was done in the grazing collision with nonzero 



impact parameter b 17|, ll8|], the upper bound on the radiation efficiency is very weak for 



large b. Hence, the direct calculation of gravitational radiation in the grazing collision is 
more important than the head-on collision case. Also, it is required to calculate the Kerr 
parameter of the resultant black hole. In order to solve these problems, the off-axis collision 
of the Bowen-York black holes could be an approximation and provide several implications. 
We plan to study this process using the close-slow approximation as a first step. 
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APPENDIX A: OTHER SOLUTIONS OF EXTRINSIC CURVATURE 

We briefly comment on other solutions of K^ v for the momentum constraint (j7]), which 
were not introduced in Sec. II. Since the equations for and x are pure Laplace equations, 
we can make infinite number of solutions for the extrinsic curvature. Here, we present 
solution of low-multipole momenta. 

To make a spinning black hole, we write 

(n - \)irGJ lw n v 
B » = ^KR- ' x = (A1) 

where J^ v is an anti-symmetric tensor. Then, K^ v is written as 

Kpu = - ^^nll ( J w nPn » + J» P n p n^) . (A2) 

J^v denotes an angular momentum tensor expressed by 

1 r 

J {w = q ^7 / {x^Kup x u Kp j p)n p dS. (A3) 

JR->oo 

In the three-dimensional case, we can define the angular momentum vector by the formula 

J, = \e,u P r p - (A4) 

The conformal factor & can be determined by solving the Hamiltonian constraint in the 
same procedure as in Sec. III. 

It is possible to derive other solutions, e.g., by setting 

5, = ^, X = 0, (A5) 
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or 



B„ 



C, 



R n+1 



n v n p - 



n + 1) 



(A6) 



where C^ u and C^ up are symmetric and anti-symmetric tensors, respectively. However the 



physical meaning of these solutions is unclear. 



APPENDIX B: INITIAL CONDITION OF THE MASTER VARIABLE 



In this appendix, we explain how to set the initial condition of the master variable $. 
We begin by studying the initial condition of the following perturbative metric 

( fH S Fi§ h E e fcnSwL \ 

sym. f~ 1 H 2 E> hi§g 
sym. sym. 2r 2 (H L ^ ee S + H T §> e 
\ sym. sym. sym. 2r 2 (if L 70 t( ^§ + H T §<p r fa) J 

where § denotes the hyper-spherical harmonics 

S = S l r ] =Kl n] Cl in - 1)/2] (cos6), 



5h 



2r 2 {H Ll<tH ^ + H T % 



47r (n+l)/2 F ( n + £ _ -g 



(n + 2t - 1)T(£ + l)r((n - l)/2)r(n - 1) 
on the unit sphere of which metric is 



-1/2 



^ ij dz l dz j = d6 2 + sin 2 OdVL 2 ^. 



Sj and SL- are given by 



As 



AS, 



(Bl) 



(B2) 
(B3) 

(B4) 
(B5) 

where A denotes the covariant derivative on the unit sphere and is defined in Eq. (JED). 



Since the initial metric is the Bri 
the ones that we have derived in 



1-Lindquist one, the values of H 2 , Hl, Ht, hi are same as 
3: 



Ho = 2Hr 



\/{n-l)R n - 1 



l + l/AR! 1 - 1 \RJ v "" 2 
h 1 = H T = 0. 

Here our purpose is to compute the time derivative of these variables. 



(B7) 
(B8) 
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The nonzero components of shown in Eqs. (|40p - (j42l) are rewritten as 



K 



RR 



(z P/M ) 



(n + 2' 



S 



+ n- 



R -{n+l) 



{z P/M 



U + 2S? , 2(n+l) g^- , 



1 id 



^-(n-l)_ 



(B9) 



(BIO) 



Thus Kn U is composed of i = and 2 modes. Since £ = mode is absorbed in a coordinate 



transformation, we omit it hereafter. The extrinsic curvature K 

in the (r, 9, 0*) coordinate is 



[where !Po is defined in Eq 

K rr = -(z P/M 



(n + 2) 



j-l^-2(n+l)/(n-l)^-(„+l) 



n + 2 # 2(n + l)4 n : 

\n\ + 



is written as 



K„ 



1 



n-1 K W 



- IV Mf 



(Bll) 



(B12) 



(B13) 



^ 2a 

in terms of the lapse function a = \fj and the shift vector /3 M . Although can be freely 
chosen provided that 0(f3^) = O(zq), we impose (3^ = for simplicity. This brings us 



#i = fro 



and fr 



-2^/fK^. Using Eqs. (jBTlil and rfB12l . we have the explicit 



formulas of h rr and and they are compared with Eq. (1B1I) : 

fri = 0, 

# 2 = 2(n + 2)?7(r), 
n + 2 



where 



Hf = 
rj(r) - 



n 

2(n + l 



r?(r), 
7](r) 



n — 1 

(ZqP/Mq) y/J 



(B14) 

(B15) 
(B16) 

(B17) 
(B18) 



Now we are in a position to derive the equation for $ and Setting 



ah 



fH #i 
#i #2// 



r/ a = (fro,frl) 



(B19) 
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the gauge invariant quantities derived in 12J are 



F = H L + (l/n)H T + (l/r)D a rX a , 



F a b — fab + D a Xb + D X a 



(B20) 
(B21) 



where 

X a = T - (/« + ^D a H T ) , (B22) 

and D a denotes the covariant derivative with respect to the metric g a bdy a dy b = —fdt 2 + 
f~ x dr 2 . In our case, 



and 



X t — JpH-T, X r — 0, 



F = X (r)/2. 



(B23) 



(B24) 



F is related to $ and $ r as shown in [12[| and we find the same equation for $ as that in our 

n ' 

previous paper [6|. Hence, the initial value of $ is unchanged. Calculating F rt , we obtain 



2(n + l)r 2 



rt 



r f 



ij(r) + 2r) >r (r) 



n — 1 k 2 

where we have substituted Eq. (1B17j) . Since F rt is related to $ and $ jr as 

Pz 



(B25) 



F rt = r 



l-n/2 



AHf 



we find the equation 



Pz 



2n 
n — 1 



,n/ 2 -l 



7y(r) 



(B26) 



(B27) 



where Pz is a function of r given in [12| and we used Eq. ( 1B18I) . The solution is given by 
Eq. (US}. 
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